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COMPUTING FAITHFUL REPRESENTATIONS 
FOR NILPOTENT LIE ALGEBRAS 

DIETRICH BURDE, BETTINA EICK, AND WILLEM DE GRAAF 

Abstract. We describe three methods to determine a faithful representation of 
small dimension for a finite-dimensional nilpotent Lie algebra over an arbitrary 
field. We apply our methods in finding bounds for the smallest dimension fj,{g) of 
a faithful g-module for some nilpotent Lie algebras g. In particular, we describe 
an infinite family of filiform nilpotent Lie algebras f„ of dimension n over Q and 
conjecture that /i(fn) > n+1. Experiments with our algorithms suggest that /Lt(fn) 
is polynomial in n. 



1. Introduction 



The Ado-Iwasawa theorem asserts that every finite-dimensional Lie algebra over 
an arbitrary field has a faithful finite-dimensional representation. A constructive 
proof for this theorem in characteristic has been given in [5] . It has been imple- 

■^' I mented as an algorithm in the computer algebra systems GAP and Magma. 

^ ' Here we consider a variation on this theme: we introduce three algorithms for com- 

puting a faithful finite-dimensional representation for a finite-dimensional nilpotent 
Lie algebra over an arbitrary field. They take as input a nilpotent Lie algebra q 



00 ! given by a structure constants table and can be briefiy summarized as follows. 

• Our first algorithm uses an action of g on a quotient of the universal enveloping 
algebra of g. It is based on ideas in |2j. 

/\ • • The second algorithm constructs a finite-dimensional faithful submodule of the 

j^ ■ dual of the universal enveloping algebra of g. The resulting module is minimal 

in the sense that it has no faithful quotients or submodules. 

• The third algorithm uses a randomised method to try to construct a faithful 
representation in dimension dim(g) + 1. It uses induction on a central series 
of g, in each step extending the representation by means of a cohomological 
construction. 

The methods described here are practical. In particular, all three methods are 
usually more efficient than the general algorithm given in [5j. We include a report 
on various applications of our GAP implementation of our algorithms below. 
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A central aim in the design of our algorithms was the idea to try to construct a 
faithful representation of possibly small dimension of a given nilpotent Lie algebra. 
In particular, our three algorithms usually determine representations which are a lot 
smaller than the representation obtained by the algorithm from [5j . Our algorithms 
can therefore be used to determine good upper bounds for the minimal possible 
dimension /i(0) of a faithful module of the nilpotent Lie algebra g. 

The interest in the invariant fi{Q) is motivated, among other things, by problems 
from geometry and topology. For example, Milnor and Auslander studied general- 
izations of crystallographic groups; these are related to /^(g) as follows: let F be a 
finitely-generated, torsion-free nilpotent group of rank n with real Malcev comple- 
tion Gr- Then, if F is the fundamental group of a compact complete afiinely-flat 
manifold, it follows that /i(0r) < n + 1 for the Lie algebra gr of Gr- Another moti- 
vation for studying /i(g) is based on the result that the Lie algebra g of a Lie group 
G admitting a left-invariant affine structure satisfies /^(g) < n + 1. In the context 
of these two results Milnor conjectured that any solvable Lie group should admit a 
left-invariant affine structure. The algebraic formulation of this conjecture implies 
that fi{g) < dim(g) + 1 should hold for the Lie algebra g of G. 

There are counter-examples known to Milnor's conjecture. Indeed, there are in- 
finitely many filiform nilpotent Lie algebras of dimension 10 which do not have any 
faithful module of dimension n + 1. We refer to [5] for details and background. 

Here we describe an infinite family of Lie algebras f„ of dimension n for n > 13 
and we use our algorithms to study the invariant /i(fn) for these Lie algebras. We 
conjecture that these Lie algebras do not have a faithful representation of dimension 
n + 1. But our experiments suggest that yu(fn) is poljTiomial in n for these Lie 
algebras. 



2. Using quotients of the universal enveloping algebra 

Let g be a finite-dimensional nilpotent Lie algebra over an arbitrary field. By g™, 
m > 1, we denote the terms of the lower central series of g. If xi, . . . , x^ is a basis 
of g, then the formal products x"^ ■ ■ ■ x'^'^ with Oj G N form a basis of the universal 
enveloping algebra f/(g). 

We define the weight wgt{x) of an element x G g as the maximal m with x G g"^. 
The weight of a basis element of U{q) is then defined by 

d 
wgt{x"^ ■ ■ ■ x^'') = ^ Oi wgt{xi) 

i=l 

Let f/™(g) = (x° I wgt^x"') > m) the ideal in f/(g) generated by all basis elements 
of weight at least m for some m>l. The following theorem is proved in [1] , [2] • 

Theorem 2.1. // g is nilpotent of class c, then g acts faithfully on U{q)/U'^^^{q) 
by multiplication from the left. If the considered basis of q contains bases for g™ for 
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every m > 1, then the resulting representation has the dimension 



iy{d,c) = J2y^_^^jPU)^ 



j=0 

where p{j) is the number of partitions of j with p{0) = 1. 

This theorem yields a straightforward algorithm to construct a faithful module 
for Q. We consider a basis of g which contains bases of g"^ for every m. We form 
the space V spanned by all basis elements of weight at most c in U{q). An element 
a; e g acts on V by left multiplication, where we treat any element of weight at least 
c + 1 as zero. We demonstrate this algorithm in the following example. 

Example 2.2. Let g be the 2i- dimensional Heisenberg Lie algebra spanned by x, y, z 
with non-zero bracket [x,y] = z. The nilpotency class of g is 2. We form the space 
spanned by the basis elements of f/(g) of weight at most 2. These are 

l,x,y,z,x'^,xy,y^. 

Thus we obtain a 7-dimensional representation for q. It is straightforward to deter- 
mine this representation explicitly by computing actions. For example, y ■ x = yx = 
xy — z and x ■ xy = x'^y = 0. 

The modules resulting from this simple and very efficient algorithm still have 
rather large dimension. In the remainder of this section, we describe two methods 
to determine a module of smaller dimension from this given module. As a first step, 
we note that a nilpotent Lie algebra g acts faithfully on a module if and only if its 
center Z{g) acts faithfully. Thus if / is an ideal in ^(g) such that / fl Z{g) = 0, 
then L acts faithfully on U{g)/I. 

For our first method we assume that the considered basis of g additionally contains 
a basis for the center Z{g). We wish to determine an ideal / in U{g) which has 
possibly small codimension and satisfies that I H Z{g) = 0. Let B be the set of 
all basis elements of weight at least c + 1 in U{g) and initialise I = (B). We now 
iterate the following procedure: let a be one of the finitely many basis elements of 
f/(g) not contained in / and not contained in Z{g). li xa E I for all x E g, then we 
add a to i? and thus enlarge / without destroying the property / fl Z{g) = 0. This 
approach usually yields a rather small dimensional faithful representation of g. We 
demonstrate this in the following example. 

Example 2.3. We continue Example \2.B, We initialise B as the basis elements of 
weight at least c+ 1. Note that Z{g) = {z). Thus we consider for a the elements 
x,y,x^,xy,y^. Of those, the elements x"^ , xy and y"^ satisfy the condition of the 
algorithm and thus we move these into B. Now also y satisfies the condition and 
we also move y into B. Now I is an ideal of codimension 3 in U{g) and hence we 
obtain a 3- dimensional faithful representation of g. 

So we now have a procedure to construct a possibly small set of monomials, that 
yields a faithful g-module: we first carry out the procedure of Theorem 12.11 and 
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then the procedure outhned above, to make the set of monomials smaller. In the 
sequel this algorithm will be called Regular. 

A second method, that can be tried for any faithful g-module V, is to perform 
the following algorithm on V. 

(1) Compute the space S = {v E V \ x ■ v = ioi all x E g}. 

(2) Compute the space C = {x-v\veV,xE Z{q)}. 

(3) Set M = SnC and let W he a complement to M in 5*. 

(4) li W = then the algorithm stops, and the output is V. Otherwise, set 
V := V/W, and return to (1). 

We note that any subspace of S* is a g-submodule of V. Therefore the quotient 
V/W is a 0-module. Let f/ be a complement to PF in \^ such that C G U. Then 
X ■ V^ C t/ for all X G Z{q). So since Z{q) acts faithfully on V it acts faithfully on 
V/W. Hence V/W is a faithful 0-module. 

Example 2.4. We consider the module of Example \2.B. Here we get 

S = {z,x^,xy,y'^) 

C={z) 

W = {x'^,xy,y'^). 

After taking the quotient we get a module spanned by (the images of) l,x,y,z. For 
this module we can perform the algorithm again. We get S = {y,z), C = (z), 
W = (y). So we end up with a faithful module of dimension 3. 

Now the complete algorithm to construct a small-dimensional faithful g-module 
consists of first performing the algorithm Regular, followed by the quotient procedure 
described above. This algorithm will be called Quotient. 

3. Using the dual of the universal enveloping algebra 

Now g acts on the dual U{g)* by x- f{a) = f{—xa). Let zi, . . . ,Zrhe a basis of the 
center of g, which we assume to be a subset of the basis xi, . . . , x„. Let ipi ^ U{q)* 
for 1 < 2 < r be defined by ipi{zi) = 1 and '?/'i(a) = for any PBW-monomial not 
equal to Zi (note that this definition depends on the choice of basis of g). 

Let " : f/(g) — > U{q) be the antiautomorphism induced by x = —x for x G g. 
Then for a,b G f/(g), / G f/(g)* we get a ■ f{b) = f{ab). (In other words, " is the 
antipode of U{q).) Note that a = a. 

Theorem 3.1. Let V be the g-submodule of U{g)* generated by ipi, . . . ,ipr- Then V 
is a faithful finite- dimensional Q-module. Moreover, V has no faithful g-submodules, 
nor has it faithful quotients. 

Proof For A; > 1 let f/'=(g) be as in Theorem O Let ly = {/ G f/(g)* | 
/(f/^^^(g)) = 0}. Then W is finite-dimensional (since t/^^^(g) has finite codimen- 
sion), and a g-submodule of U{q)* (since U'^~^^{g) is an ideal). Now V C W; hence 
V is finite-dimensional. Let z = Yli f^i^i be an element of the center of g. Then 
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z ■ ipill) = —ipilz) = —^i- So z acts as zero if and only if all /Xj are zero. So V^ is a 
faithful module. 

Let TpQ G U{q)* be defined by ^o(l) = 1 and ipQ^a) = for PBW- monomials a 7^ 1. 
Then ^/'o = ~zi ■ ipi, so ipo & V . Let f & V, and suppose there is a PBW- monomial 
a 7^ 1 such that /(a) 7^ 0. Then a ■ /(I) = /(a) 7^ 0. The conclusion is that ipo 
spans the space of elements that are killed by g. Set M = V/ {ipo). Then M is not 
a faithful g-module. Indeed, V has a basis consisting oi a ■ ipi for 1 < i < r, and 
various PBW-monomials a. Let z lie in the center of q; then for all PBW-monomials 
b we get z ■ {a ■ ipi){b) = ipiiabz) which is zero unless a = h = 1 (note that z = —z 
also lies in the center of q). Furthermore, z ■ ipi = —fiiipo (where z = J^jf^j^j)- ^^ 
follows that the center of q acts trivially on V/ {ipo). In particular it is not a faithful 
g- module. Now, since every g-submodule of V must contain ■j/'q, it follows that V 
has no faithful quotients. 

Let M C y be a faithful g-submodule. Let hij E U{q) be such that {bijipi} is a 
basis of V. We assume that bn = 1, and that the bij G qU{q) if j > 1 (i.e., they 
have no constant term). Then Zi ■ buipi = —ipo and Zi ■ bkjipk = if /c 7^ 2 or j > 1. 
So since the center acts faithfully on M it follows that M contains elements of the 
form 

V^i = ^i + X] ^kjbkji'k, 
i>i 

l<fc<r 

for 1 < i < r. (Here Ckj are coefficients in the ground field.) Now we introduce a 
weight function on [/(g)*. For fc > set Ft = {f E [/(g)* | /(f/^'(0)) = 0}, where 
U''{q) is as in Theorem EH Then = Fq C Fi C ■ • • . We set wgt{f) = k ii f E Fk 
but / ^ Fk-i- (For example: wgt{ipo) = 1-) Let / G U{q)* have weight k, and let 
a E U{q) with wgt{a) = t; then a calculation shows that wgt{a ■ f) < k — t. Hence 
bij(fi is equal to bijipi plus a sum of functions of smaller weight. So if we order the 
bijipi according to weight, and express the bijipi on the basis hjipi we get a triangular 
system. We conclude that the bij(fi are linearly independent. Hence dimM = dim V^ 
and M = V. D 

The algorithm based on this theorem is straightforward. We illustrate it with an 
example. 

Example 3.2. Let g be the Lie algebra of Example {KM For a monomial a E U{g) 
we denote by ipa the element of U{q)* that takes the value 1 on a, and zero on all 
other monomials. We compute a basis of the submodule ofU{g)* generated by ipz- 
We have x ■ ipz{ci) = —i^zi^a) = for all monomials a. Secondly, y ■ ipz{x) = 
Tpz^—yx) = ipzi—xy + z) = 1. So we get y ■ ipz = i^x- Furthermore, z ■ ipz = —''Pi 
and Q ■ ipi = Q, X ■ ipx = —''Pi, y ■ ''Px = z ■ ipx = 0. So the result is a ?>- dimensional 
g-module. 

Remark 3.3. Since we are working in the dual of an infinite-dimensional space it is 
not immediately clear how to implement this algorithm. We proceed as follows. Let 
V be as in Theorem 13. II From the proof of Theorem 13. II it follows that f{U{QY~^^) = 
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for all / G y. In other words, for all monomials a with wgt{a) > c + 1 and all 
f & V we have f{a) = 0. It follows that we can represent an / G V by the vector 
containing the values /(a), where a runs through the monomials of weight < c. This 
enables us to perform the operations of linear algebra (testing linear dependence, 
constructing bases of subspaces and so on) with the elements of V. Furthermore, 
we can compute the action of elements of g on V. 

The disadvantage of this approach is that the number of monomials that have to 
be considered can be very large. So in the same way as in the previous section we 
try to throw some monomials away. Let A be the set of monomials relative to which 
we represent the elements of U{q)* . At the start this will be the set of monomials of 
weight < c. Let B be the set of all other monomials. So at the outset B spans a left 
ideal of U{g) and f{b) = for all / G \^ and b E B. We move elements from A to 
B, without changing this last property. Let a E A he such that a ^ Z{q) and xa is 
a linear combination of elements of B for all x E g. Then we claim that /(a) = for 
all / G V^. In order to see this we use the basis {bijipi} used in the proof of Theorem 
13.11 If j = 1 then bij'ijji{a) = 'ipi{a) = as a ^ Z{g). If j > 1 then bij G gU{g) and 
hence bija is a linear combination of elements in B. Hence bijipila) = 0. Also the 
span of B along with a continues to be a left ideal. We conclude that we can move 
a from A to B. We continue this process until we we do not find such monomials 
any more. The resulting set is usually a lot smaller than the initial one. 

We note that the procedure described in the previous remark is exactly the same 
as the second phase of the algorithm Regular (see Section [2]). So we first perform 
the algorithm Regular, and use the resulting set of monomials to represent elements 
of the dual of U{g). The resulting algorithm is called Dual. 

4. AfFINE REPRESENTATIONS AT RANDOM 

Let g be a nilpotent Lie algebra of dimension d. A homomorphism p : q ^ 
aff{K'^) C g[d^i[K) into the Lie algebra of affine transformations 

is called an affine representation of g. In this section we describe a method that tries 
to determine a faithful affine representation of g of dimension d + 1. If the method 
succeeds, then it returns such a faithul representation of dimension d+1. However, 
it may also happen that the method fails and does not return a representation. Also, 
it is worth noting that the algorithm uses random methods and hence different runs 
of the algorithm may produce different results. 

The method uses induction on a central series in g. Thus we assume by induction 
that we have given a central ideal / in g with dim(/) = 1 and a faithful affine 
representation 

p:g/I^MdiK). 

Let {ai, . . . , ttd} be a basis of g with I = (ad) and let Mj = p{ai + I) for 1 < 2 < d — 1. 
We assume that every Mj is a lower triangular matrix. Clearly, we can readily extend 
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p to an afiine representation of Q with p{ai) = Mi for 1 < i < d where we set M^ = 
(so that p(a(i) = 0). This extended representation has kernel /. 
Our aim is to extend p to a faithful afiine representation 

iP : Q ^ Md+i{K) 

such that 

Ha^) =(^' qO for 1 < ^ < d, 

for certain vectors Vi G K'^. The following lemma shows that the possible values for 
Vi can be determined using a cohomology computation. Recall that 

Z\q,K'') = {u:g^K'' linear | u{[x,y]) = p{x)u{y) - p{y)u{x)} 

is the space of 1-cocycles with values in the p(0)-module K'^. 

Lemma 4.1. ip is a representation of q if and only if Vi = 6{ai) for 1 < i < d for 
some 6 G Z^{g,K'^). 

Proof. Let 6 G Z^{q,K'^) with 6{ai) = Vi. The linearity of 5 implies that ip is 
linear. The defining condition for maps in ^^(0, K*^) implies that ■?/' is a Lie algebra 
representation. The converse follows with similar arguments. D 

Note that Z^{q, K'^) is a vector space over K and can be computed readily using 
linear algebra methods. The computation of Z^i^Q^K*^) allows to describe all afiine 
representations of q extending p. It remains to determine the faithful representation 
among these. 

Lemma 4.2. ip is faithful if and only if Vd+i 7^ 0. 

Proof. If ip is faithful, then Vd+i 7^ 0. Conversely, suppose that Vd+i 7^ 0. As p is 
faithful, it follows that ker{ip) C /. As Vd+i 7^ 0, we find that ker{ip) = 0. D 

These ideas can be combined to the following algorithm. 

(1) Choose a central series S = go > 0i > • • • > 0d > 0d+i = of ideals in g such 
that dim(gi/0j+i) = 1. 

(2) by induction, extend a faithful afiine representation from g/gj to g/gj+i: 

- Compute Z^{g/Qi+i,K'). 

- Choose a 5 e Z^{Q/gi^i, K^) with 6{ai) 7^ 0. 

- If no such 6 exists, then return fail. 
-US exists, then extend p to g/g^+i. 

If g has a faithful afiine representation of dimension d+1, then this algorithm can 
in principle find it. However, it may be that a "wrong" choice of a 5 at a certain 
step may cause the algorithm to fail at a later step. 

The algorithm is based on linear algebra only and hence is very effective. It often 
suceeds in finding a faithful representation in dimension d + 1 if it exists. 
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5. A SERIES OF FILIFORM NILPOTENT LlE ALGEBRAS 

Let Khe a field of cliaracteristic zero. In tliis section we define a filiform Lie algeba 
fn in each dimension n > 13 having interesting properties concerning Lie algebra 
cohomology, affine structures and faithful representations. In fact, we believe that 
the algebras f„ are counter examples to the conjecture of Milnor mentioned in the 
introduction, i.e., that /i(f„) > n + 2 holds. Hence it is interesting to compute the 
invariants /u(f„). 

Define an index set X„ by 

jO = {(A;, s)gNxN|2<A;< [n/2], 2A; + 1 < s < n}, 
1^ if n is odd, 

n 
2- 



^2^U{(|,n)} if n is even 



Now fix n > 13. We define a filiform Lie algebra f„ of dimension n over K as 
follows. For (/c, s) G X„ let ak^s be a set of parameters, subject to the following 
conditions: all ak s are zero, except for the following ones: 



3 

Oie,2£+i = ,n /2f-n ' -^ = 2, 3, . . . , L^J ■, 

\2l U-1 j 

_ 1 10(n-7)(n-8) 
"^•"-2"7 + 21(n-4)(n-5)' 

22105 if„-1Q 

' ^0 if n > 14, 



and 

1 _ 70(n-8) 25(n-6)(n-7)(n-8) 

«5,n - ^ - ii(^_2)(n-3)(n-4)(n-5) ^ 99 (n - 2)(n - 3)(n - 4) 

5 (n-5)(n-6) 65 (n-7)(n-8) 



+ 



66(n-2)(n-3) 1386 (n - 4)(n - 5)' 



Let (ei, . . . , e„) be a basis of f„ and define the Lie brackets as follows: 
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[ei,ei] = Cj+i, i = 2,...,n-l 



[Cj, CjJ 



J-t-i-l 



X^f X] (-l)M p )«i+^,r-i+i+2^+i)er, 2<i<j<n. 



In order to show that this defines a Lie bracket we need the following lemma 
which follows from the Pfaff-Saalschiitz sum formula: 

Lemma 5.1. We have the following identities for all n > 13; 



> , (-1) n 2 , "^'2^+1 = 
i=3 ^ ^ 


(n-7)(n-8) 




(n-4)(n-5)' 




1 "~^ 
£=5 ^ ^ 


1 12(n-8) 




70 ' (n-2)(n-3)(n-4)(n- 


-5) 


1 "~^ 


{n-5){n-6) 
(n-2)(n-3)- 





Proposition 5.2. The Jacobi identity is satisfied, so that f„ is a Lie algebra for any 
n> 13. 

Proof. Let n > 14 and choose the parameters a^^s as follows. Consider ak,2k+i, 
k = 3, . . . , |_^^J and a4^n-2, ct5,n as free variables. Let the remaining parameters be 
zero, except for 02,5 = 1, ct3,7 7^ and q;3,„_4 = 1. The Jacobi identity is equivalent 
to a system of polynomial equations in the free parameters. First we obtain the 
equation a^ji^lOasj — 02,5) = 0, so that 0:3,7 = j^. More generally we see that 

(£ - 1) ■ a,,2£+i = (4£ + 2) ■ a£+i,2£+3, ^ = 2,3,..., [^J . 
This implies the given explicit formula for all 0^,2^+1- Secondly we obtain 



L^J 



2 
0:4,9 

OiAn-2 






'^3,7 0.1 £^3 
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Otb,n 




H \ 7> «4,n-2 13a4,9+ y] (-1)'^( f, o )««.2£+l 

"4,9 + "3,7 - 2a2,5 I \ f:^ \ i-2 J 



i=3 



This amounts to the given formulas in the definition of f„, if we substitute the 
identities from Lemma 15.11 Conversely this also shows that the Jacobi identity is 
satisfied if the free parameters are given in this way. 

For n = 13 there is one difference. The parameter ^4 „ coincides with the parameter 
0:4 13, which is given by 

_ «3,9(-«5,13 + 604,11 - 503,9) 

'^4,13 — -7^ , 

"3,7 + 2^2,5 

and cannot be chosen to be zero. For ra > 14 the choice a^^n = is consistent with 
the Jacobi identity. D 

Example 5.3. The parameters for fi3 are given as follows: 

_ 1 _ 1 _ 1 1 

0^2,5 — I5 Oisj — TTT) ct4,9 — t;^: Oi5,ii — 77^5 Oie,13 



10' '"" 70' '""^^ 420' '"''" 2310' 



43 22105 313 

«3,9 = 1, "4,11 = 7737, «4,13 = 7777775 «5,13 



126 ' 15246 ' 3388 

The algebras f„ belong to the family of filiform Lie algebras 21^(7^) defined in [3]. 
Let us recall the following definition. 

Definition 5.4. Let g be a filiform nilpotent Lie algebra of dimension n. A 2- 
cocycle u G Z^{q, K) is called affine, if cu: g A g — > i^ does not vanish on 3(g) A g. 
A class [uj] G H^{q, K) is called affine if every representative is affine. 

The cohomology class [uo] G H^{q, K) of an affine 2-cocycle uo is affine and nonzero. 
If a filiform Lie algebra g of dimension n > 6 has second Betti number 62(5) = 2, 
then there exists no affine cohomology class. 

We have shown in [5] that a filiform Lie algebra g which has an affine cohomology 
class, admits a central extension 

0^a^[)^g-^0 

with some Lie algebra () and i(a) = ^(P)), and has an affine structure. In particular, 
such a Lie algebra has a faithful representation of dimension n + 1. 
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We can conclude from the results in [3] that the Lie algebras f„ do not have an affine 
structure arising this way. 

Proposition 5.5. The algebras f„, n > 13 have second Betti number &2(fn) = 2. 
Hence there exists no affine cohomology class [u] G H^{q,K). In particular there is 
no central Lie algebra extension as above. 

For Lie algebras in 2t,^(i^) the second Betti number is 3 or 2, depending on whether 
a certain polynomial identity a3,n-4 = -P in the free parameters does hold or does 
not hold. For f„ we have chosen the parameters in such a way that P = and 
Q^3,n-4 = 1- This implies that 62(fn) = 2. 

It follows that a very natural way to obtain a faithful representation of dimension 
n + 1 does not work. In fact, we believe that there is no such representation at all 
for these algebras: 

Conjecture 5.6. The Lie algebras f„, n > 13 do not have any faithful representation 
of dimension n + 1, i.e., /i(fn) > n + 2. 

For n = 13 a very complicated analysis of possible faithful representations seems 
to confirm this conjecture. In general our methods are not sufficient to prove this 
for all n > 14. Even more difficult of course is the determination of /u(f„). 

6. Practical experiences 

We implemented all the algorithms described above in the computer algebra sys- 
tem GAP. In this section we report on the application of these implementations to 
various examples. From Section [2] we have the algorithms Regular and Quotient. 
From Section [3] we have the algorithm Dual. Finally the algorithm of Section H] is 
called Affine. 

In all our experiments Quotient and Dual returned faithful representations of the 
same dimension (with Dual being slightly faster). This is illustrated in Table O 
We believe that there must be an intrinsic reason for this to happen, such as one 
module being the dual of the other. But we have no proof of that. We only exhibit 
the results of Dual in Tables [1] and [2], noting that the results for Quotient are similar 
in all cases. 

All computations were done on a 2GIIz processor with 1GB of memory for GAP. 

6.1. Upper triangular matrix Lie algebras. The upper triangular matrices in 
M„(F) form a nilpotent Lie algebra t/„(F) with n — 1 generators and class n — 1. We 
applied our algorithms to some Lie algebras of this type. The results are recorded 
in Table 1. 

Table 1 exhibits that the underlying field does not have much impact on the 
runtime or the result. The larger the dimension of the considered Lie algebra is, 
the more superior is Affine. It yields small dimensional representations and is the 
fastest of all methods. 
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n 


F 


dim(t/„(F)) 


Regular 


Dual 


Affine 








time 


dim 


time 


dim 


time 


dim 


4 


F2 


6 


0.0 


7 


0.1 


5 


0.0 


7 


5 


F2 


10 


0.25 


15 


0.3 


11 


0.3 


11 


6 


F2 


15 


3.4 


35 


3.6 


17 


3.5 


16 


7 


F2 


21 


65 


79 


66 


35 


45 


22 


4 


F3 


6 


0.0 


7 


0.0 


5 


0.0 


7 


5 


F3 


10 


0.2 


15 


0.3 


11 


0.3 


11 


6 


F3 


15 


3.4 


35 


3.6 


17 


3.7 


16 


7 


F3 


21 


65 


79 


67 


35 


46 


22 


4 


Q 


6 


0.0 


7 


0.0 


5 


0.0 


7 


5 


Q 


10 


0.2 


15 


0.3 


11 


0.3 


11 


6 


Q 


15 


3.0 


35 


3.2 


17 


3.6 


16 


7 


Q 


21 


66 


79 


67 


35 


45 


22 



Table 1. Running times (in seconds) for f/„(F). 



6.2. Free nilpotent Lie algebras. Next we consider the free nilpotent Lie algebras 
with n generators of class c over the field F, denoted N^ c{^)- 



n 


c 


F 


dim(iV„,e(F)) 


Regular 


Dual 


Affine 










time 


dim 


time 


dim 


time 


dim 


2 


5 


Q 


14 


0.2 


20 


0.3 


20 


0.5 


15 


2 


6 


Q 


23 


0.9 


34 


1.3 


34 


8.4 


24 


2 


7 


Q 


41 


3.2 


65 


4.8 


65 


© 


© 


2 


8 


Q 


71 


14 


117 


21 


117 


© 


© 


3 


4 


Q 


32 


0.8 


41 


1.7 


41 


54 


33 


3 


5 


Q 


80 


11.5 


113 


17.5 


113 


© 


© 


4 


3 


Q 


30 


0.9 


36 


1.3 


36 


37 


31 


4 


4 


Q 


90 


13 


113 


19.7 


113 


© 


© 



Table 2. Running times (in seconds) for N^ 



Table 2 displays the time in seconds for the three algorithms, with input A^„,c- The 
© in the last two columns indicates that the algorithm Affine did not succeed, either 
because it made the "wrong" choice at some stage, or due to Memory problems: for 
its cohomology computation it has to solve a system of linear equations which is of 
the size 0(dim(g)^) and this can be time-and space consuming. 

6.3. The Lie algebras f„. Finally, we consider the Lie algebras f„ of the previous 
section. The results of that are contained in Table O 
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n 


Regular 


Quotient 


Dual 


Affine 




time 


dim 


time 


dim 


time 


dim 




13 


8.6 


85 


14 


43 


12.3 


43 


© 


14 


17 


105 


28 


53 


24.7 


53 


© 


15 


33 


145 


63 


64 


50 


64 


© 


16 


64 


185 


125 


i i 


102 


77 


© 


17 


123 


256 


323 


94 


218 


94 


© 


18 


234 


316 


731 


111 


461 


111 


© 


19 


487 


433 


1844 


134 


1162 


134 


© 


20 


920 


538 


4009 


158 


3039 


158 


© 



Table 3. Running time (in seconds) for the Lie algebras f„. 



Table 3 displays the time in seconds for the algorithms Quotient and Dual, with 
input f„. The © in the last column indicates that the algorithm Affine did not 
succeed. In this case, this was due to the fact that Affine did not find any possible 
faithful representation of dimension n + 1. Of course, if our conjecture on f„ holds, 
then it cannot succeed. 

Note that the dimensionals of the determined modules for f„ are significantly 
larger than n + 1. However, they do not seem to grow very fast. Some naive tests 
with least squares fits seem to suggest that the dimensions grow quadratically or 
cubically. 

6.4. Some comments. From the above tables we conclude that if Affine succeeds, 
then it usually finds a module of significantly smaller dimension than Regular, Quo- 
tient or Dual. This supports the suggested strategy to try this algorithm first. 

References 

[1] G. Birkhoff: Representability of Lie algebras and Lie groups by matrices. Annals of Mathe- 
matics 38 (1937), 526-532. 

[2] D. Burde: A refinement of Ado's Theorem. Archiv Math. 70 (1998), 118-127. 

[3] D. Burde; Affine cohomology classes for filiform Lie algebras. Contemporary Mathematics 
262 (2000), 159-170. 

[4] The GAP Group: GAP - Groups, Algorithms, Programming - a System for Computational 
Discrete Algebra Version ^.^.10 (2007), http://www.gap-system.org. 

[5] W. A. de Graaf: Constructing faithful matrix representations of Lie algebras. Proceedings 
of the 1997 International Symposium on Symbolic and Algebraic Computation, ISSAC'97, 
(1997), 54-59. ACM Press New York. 



Fakultat fur Mathematik, Universitat Wien, Nordbergstr. 15, 1090 WiEN, Aus- 
tria 

E-mail address: dietrich.burdeSunivie .ac . at 



14 D. BURDE, B. EICK, AND W. DE GRAAF 



Institut Computational Mathematics, TU Braunschweig, Pockelsstrae 14, 31806 
Braunschweig, Germany 

E-mail address: beick@tu-bs.de 

Dipartimento di Matematica, Via Sommarive 14, 1-38050 Povo (Trento), Italy 
E-mail address: degraaf@scieiice.unitn.it 



